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Abstract 

The electron and phonon temperature distribution function are 
calculated in semiconductors. We solved the coupled one-dimensional 
heat-diffussion equations in the linear approximation in which the 
physical parameters on the sample are independent of the tempera- 
ture. We also consider the heat flux at the surface of the semiconductor 
as a boundary condition for each electron and phonon systems instead 
of using a fixed temperature. From this, we obtain an expression for 
electron and phonon temperature respectively. The characterization 
of the thermal waves properties is duscussed and some practical proce- 
dures for this purpose provide us information about the electron and 
phonon thermal parameters. 

1 Introduction 

Thermal wave physics is becoming a valuable tool in the study of mate- 
rial parameters as well as in the semiconductor industry for characterizing 
process in the manufacturing of electronic devices. 1-2 These waves are cre- 
ated whenever there is a periodic heat generation in a medium. The most 
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common mechanism for producing thermal waves is the absorption of an in- 
tensity modulated light beam by a sample. The photothermal heating of a 
sample leads to thermal ant stress-induced changes in the physical properties 
of both the sample and the surrounding madia. Based on this observation, 
several alternative detection techniques have been developed for monitoring 
the photothermally induced changes of given physical properties either of 
the sample or of its surrounding medium. Some of the different contact-type 
detection techniques, such as the conventional gas-microphone photoacous- 
tic detection, 1 photopyroelectric detection, 3 or remote semsing techniques as 
photothermal reflection 4 etc., have been reviewed in Ref. 5. 

Being a photothermal techniques, the detected signal is strongly depen- 
dent upon how the heat diffuses through the sample allows us to perform 
both thermal characterization of the sample (i.e., measurements of its ther- 
mal properties, such as thermal diffusivity and thermal conductivity), and 
carrier-transport properties. 6 In the case of semiconductors the photoacous- 
tical signal can provided us with additional information regarding the heat- 
transport properties and electron-phonon energy relaxation by the interact- 
ing electron and phonon systems, a fact which has been recognized since the 
earlier theory of heat conduction in semiconductors and hot electrons 6 8 
Qualitatively, this may be understood as follows: in studies of the behavior 
of semiconductors heating by electromagnetic waves it is usual to assume 
that the phonon system remains in equilibrium, i.e., that the phonon tem- 
perature T p is equal to the temperature T of the ambient medium. 8 However, 
nonequilibrium carriers in the bulk of the sample can interact and transfer 
energy to the phonon system producing an essential increase in the average 
energy of phonons (energy nonequilibrium), which can be most conveniently 
described in terms of heating of the phonon gas (increase of its temperature 



Let v pe (q) be the characteristic phonon relaxation frequency owing to 
interaction with electrons. Gurevich at al. 9 showed that the value of v ve de- 
creases rapidly for phonon wavevector q > 2p, tending ro zero, where p is 
the average electron momentum, namely p ~ \/2mT e for the nondegenerated 
electron gas, and p ~ pp (Fermi momentum) for the case of degeneracy. 
This mean that electrons are scattered by phonones having momentum val- 
ues smaller than the average electron momentum (q < 2p). Those phonones 
are knowen as long- wave (LW), in contrast to short-wave (SW) phonones, 
whose momentum satisfy q > 2p. The degree of nonequilibrium of the LW 
(and hence SW) phonones is determinated by the relationship of the phonon- 
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phonon collision frequency v vv and v pe : 
i) Strong phonon-phonon interaction 

This case satisfies the following inequality v pp ^> is pe . This situation means 
that the phonon-phonon collisions are more frequent than phonon-electron 
collisions and more efficient in terms of energy relaxation than energy transfer 
from the phonon to the electron system. As a result, the Planck distribution 
function describes the phonon system with its own temperature T p differ- 
ent from T e . Since the phonon-phonon collision frequency depends on the 
phonon temperature, the inequality u pp ^> v pe holds for phonon temperature 
T p > 50°K for n-type Ge and GaAs semiconductors with the charge density 
concentration n oc 10 14 cm 3 . 9 

ii) Strong phonon-electron interaction 
Now, consider the situation represented by the inequality v pe ^> v pp . This 
is the case when phonon-phonon collisions alone cannot bring the phonon 
system to an internal equilibrium, and the description becomes more compli- 
cated. If T p is the characteristic phonon temperature, than the momentum 
T p / s (s is the sound velosity) sets a limitto the phase spacevolume occupied 
by phonons (with q > T p s the number of phonons is exponentially low). The 
number of LW and SW phonons depend on the relation between 2p and T p j s. 
If T p / s ^> 2p is true i.e., LW phonones occupy a much smaller phase volume 
that SW phonones, then the subsystem of SW phonones has enough time to 
redistribute the energy recived from LW phonones between its constituent 
quasiparticles. As a result, the distribution function of SW phonones becones 
Plankian, with some temperature T p w . Then the electron- LW phonon inter- 
raction relax their energy more efficiently than the phonon subsystem and 
the LW phonons emitted by electrons of temperature T e are characterized by 
the same temperature T e = T^ w . In this situation we also have two different 
subsystems; one corresponds to th SW phonons with Temperature T^ w and 
the other one corresponds to electrons and LW phonons with a characteris- 
tic temperature T e = T p LW ^ T p w . If T p /s < 2p holds, then all phonons 
interact efficiently with electrons, and hence the phonon system cannot be 
subdivided into LW and SW phonons. In this case, the phonons emitted 
by electrons of temperature T e are characterized by the same temperature 
T e = T p and all the phonons show v ve ^> v vv for T p < 50°K. 

At this point its important to mention that thermal waves phenomenon 
is present in the same only for modulated frequency of the incident light 
iv (chopper frequency) of the same order as the frequency of the relaxation 
energy between the quasiparticle systems z/ e . 8 In the limit uj v e the system 
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cannot respond to this external perturbation; therefore the dynamic part of 
the heat flux is negligible as compared with the static part and the transfered 
heat is only static. For u <C v £ , the variation of the electron and phonon 
temperature oscillates with the same modulated frequency of the incident 
light. 10 For the case when u pp u pe then u £ ~ u pe otherwise [y pp <C u pe ) 
v £ ~ v P p- 9 Typical values of v e are in the range of 10 8 -10 10 sec -1 . 6 

Recently, 10 there has been some interest in studying thermal characteri- 
zation of layered systems using the photoacoustic effect. In particular, from 
the theoretical point of view, thermal diffusion of one and two layer system 
has been investigated using the heat-diffusion equation in the approxima- 
tion when both the electron and phonon temperature distribution are equal. 
This approximation is only valid in the limit of infinite electron-phonon en- 
ergy interaction in which the size of the sample is grater than the cooling 
lenght. 7 The effective diffusivity and thermal conductivity of the two layer 
were obtained assuming continuity of heat flux, instead of the temperature 
distribution, at each interface and taking into account the physical solution 
of the dynamical part of the temperature fluctuation in each layer. 

In this work we further extend the earlier model for thin film semicon- 
ductors by taking into account the so far neglected important features the 
behavior of electtrons and phonons under time varying excitation namely, the 
electron and phonon temperature distribution function. In Part II the elec- 
tron and phonon temperature distribution for semiconductors is calculated 
by solving the coupled heat-diffusion equations, for both of quasiparticles 
systems. In Section III a discussion of our results is presented, and a com- 
parison to the predictions of the existing theories is also made. Finally, in 
Section IV we present our conclusions. 

2 Formulation 

It is well known the heat transport in solids is carried out by various quasi- 
particles (electrons, holes, phonons, magnons, plasmons, etc.). Frequently 
the interactions between these quasiparticles are such that each of these sub- 
systems can have its own temperature and the physical conditions at the 
boundary of the sample. We restrict our analysis to the case of monopolar 
semiconductors under the conditions of strong phonon-phonon interaction 
asdiscussed previously. 11 Therefore steady state heat conduction can be de- 
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scribed by the following system of equations 9 . 



div Q e = -P ep (T e - T p ), div Q p = P pe (T e - T p ). (1) 

The term P ep (T e — T p ) describes the transfer of heat between electrons and 
phonons. Here P ep is a parameter proportional to the energy frequency be- 
tween electron and phonon systems (P ep = P pe = P ~ nu £ ) and the heat 
flux of electron Q e and phonon Q p subsystems are described by the usual 
relationships: 

Q e = -x e gradT e , Q p = -x p gradT p (2) 

where x e (x p ) is the electron (phonon) thermal conductivity. So far, Eqs. (1) 
include only the static contribution of the heat transport, i.e. the heat flux is 
independent of time. However, in the photothermal experiments, the incident 
radiation is modulated in time by the chooper, and in this case it is necessary 
to consider the dynamic contribution to the heat transport in the electron 
and phonon systems. Let us assume the one dimensional model for the heat 
flux. Assuming that the sample is optically opaque to the incident light (i.e., 
all the incident light is absorbed at the surface), the electron and phonon 
temperature distribution function are the solutions of 

d 2 TJz,t) . . m . 1 dT e (z,t) 

d 2 T p (z,t) ,, rm/ . rT1/ 1 dTJz,t) 

where k 2 p = and the diffusivity for each system is given as a etP = 

x e ,p/(pc) ep and p e (p p ), c e (c p ) is the electron (phonon) density and specific 
heat respectively. 

The temperature fluctuation T e>p (z, t) should be supplemented by bound- 
ary conditions at the surface of the semiconductor (z — 0). In the photother- 
mal experiment, the most common mechanism to produce thermal waves is 
the absorption by the sample of an intensity modulated light beam with fre- 
quency modulation ou < v e . It is clear that when the intensity of the radiation 
is fixed, the light-into-heat convension at the surface of the sample can be 
written in general as 

Qe, P (^ t )\ z=0 = Qe, P + ^Qe, P e lUJt (4) 
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where Q e>p is proportional to the intensity of high frequency light and the 
other term represents the modulation of this light. 

The general solution of the coupled heat-diffusion equation for the elec- 
tron and phonon system can be written as 

k 2 

T etP = A + Bz± -^fCe~ kz + e ejP (z)e iut (5) 

where k 2 = k 2 + k 2 represents the inverse of cooling lenght 7 and 9 e>p {z) 
satisfies a similar set of Eqs. (3) but instead of the term (l/a e)P )dT e ^ p /dt we 
substitute (iu /a ejP )T ejP , and the solution is given by 

9 e {z) = Fe-° z (6) 



9 p (z) = Ge 



-<JZ 



Here F and G are related each other by the following relationship 



G = (7) 
and the values of a are given by 

*?,2 = \ {-I + °l) ± \ [(^ " -If + ^ e k 2 p ] V2 (8) 

Eq. (8) represents the condition for non-trivial solutions of the coupled elec- 
tron-phonon differential equations and o\ are given by 

<v = — + k lr ( 9 ) 



a 



e,p 



It is worth mentioning that the increasing exponential term exp kz and 
expai^z which are solution of the static and the dynamical part of heat- 
diffusion equations respectively have not been considered because they do 
not represent a physical solution (heat flux cannot be reflected) . 

Using the boundary conditions at the surface of the sample given by 
Eq. (4) the electron and phonon temperatures are given by 

k 2 

T e = A + Bz+ -^Ce- kz + e iuJt \F ie ~ aiz + F 2 e~ CT2Z ] , 
k J 



6 



k 2 

T P = A + Bz+ -^Ce~ kz + e luJt [de^ 12 + G 2 e~ a2Z ] , (10) 
where the constants B, C, F 12 , Gi,2 can be written as 

B = "K^ + ^) (lla) 



(lib) 



Fl = 



AQp 7 2 i ^Qe / 2 2\ 



;iic) 



F 2 = 



0-2(0-2 - ^i 2 ) 



^Qp 7 2 1 / 2 2\ 

k e H (o"! - a e ) 



X„ 



(lid) 



Gi >2 = -^5-^^1,2 (He) 

Here A is a constant which cannot be determinated from these boundary 
conditions and it is not important in obtaining the physical results. 

Once we know the electron and phonon temperature distributions in the 
sample, we can calculated the responce of the surrounding medium due to 
thephotothermal heating of the sample using one of the several alternative 
detection techniques mentioned before. 



3 Special Cases 

We now turn to a discussion on the results obtained so far and compare 
them with previous theories on thermal waves. We shall first consider a 
nondegenerate semiconductors. In this case, the typical ratio of the heat 
conductivity of electrons to phonon satisfies x e / ' k v ~ 10~ 3 then k e ^> k p . 
Under these circunstances, after simplying the expression for T e , T p and a 12 
can be written as 

T e (z, t) = A + Bz + Ce~ kz + e iwt (F ie ~ aiz + F 2 e- a2Z ) (12a) 
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k 2 

T p (z, t)=A + Bz- t| Ce- fe2 + G 2 e iojt ~ a2Z (12b) 



where G\ ~ 0, a"i = er e and G2 



1 - f (i - i)] ^ ^ = W«p- K is 
important to note that the contribution of the electron diffusivity dissapears 
in the phase of the dynamical part of the phonon temperature. 

Because the main source of the photothermal signal arieses from the pe- 
riodic heat flow from the semiconductor, the periodic diffusion process pro- 
duces a periodic temperature variation in the semiconductor given by the 
sinusoidal (AC) component of Eqs. (12). 

Information about electron and phonon parameters can be obtained from 
photothermal experiments depending upon the decay lenght Si ;2 = 3FJcr lj2 of 
the AC component of the electron and phonon temperature at z = d (where 
d is the thickness of the sample). 

As can be seen in Fig. la, if |Fi| ^> IG2I > \F 2 \ and si <C s 2 the time 
dependent component of the phonon temperature in the semiconductor at- 
tenuates rapidly to zero with increasing distance from the surface of the solid 
as compared with component of the electron temperature, and the electron 
thermal wave carried out all the informationabout electron and phonon pa- 
rameters through the coefficient F\. However if besides that <C 
electron thermal wave provide us only information about thermal parame- 
ters of the electron system. Now if 3> information about the phonon 
system is only obtained from the photothermal experiments. 

On the other hand if ^> \G 2 \ > I-P2I and s\ ^> s 2 , in this case the 
time dependent component of the electron temperature is fully damped out 
as shown in Fig. lb. Information about the electron or phonon parameters 
depend on the relation between d and s^ 1 and how the incident energy is 
distributed in the electron and phonon systems: 

(1) d < s' 1 and > Eqs. (11c) and (12a) tell us that thermal wave 
asociated with the electron system only give information about the electron 
thermal parameters. 

(2) d > s^ 1 and <C ^2e, in this case only a phonon thermal wave can 
be detected at z = d and the photothermal experiments give the thermal 
phonon parameters of the sample. 

Similar behavoir of photothermal waves in the semiconductors can be 
obtained if \G 2 \ ^> |i*2|, \Fi\ and s 2 <C S\ and information about phonon 
thermal parameters is obtained from the photothermal experiments if in ad- 
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ditionthe inequality > ^Q- holds. If \G 2 \ > |F 2 |,|Fi| and s 2 > s u 

>Cp >f e 11 I I I I 

hence the photothermal experiments only detect the thermal parameters of 
the phonon system when d <C s^ 1 and ^> ^— or thermal parameters of 

the electron system if d ^> s? 1 and 3> as show in Figs. 2a and 2b. 

In the case k ejP — > (electron-phonon interaction vanishes), it follows 
from Eqs. (11) that the electron and phonon temperatures reduce to 

1 /2 

Tep = Ao - Q^z + ( (1 - (13) 

X C)P 2x 6iP \ tO ) 

whete the constant A has been chosen such it exactly cancels out the divirgent 
term from the Eq. (lib), i.e. for P — > we have non-interacting systems of 
quasiparticles and heat transport is carried out by the electrons and phonon 
independently through the semiconductor. 

We now consider, the situation represented by the strong coupling be- 
tween electrons and phonons i.e. P — > oo. In this case the electron-phonon 
energy interaction is very efficient in terms of energy relaxation than energy 
transfer in the electron or phonon system. Here the solutions of Eq. (8) are 
given by 

= * + + * 



and G 2 ~ F 2 , B = — - ^ (Q e + Q p ) therefore, in this limit both tempera- 
tures are the same 



T = T e = T p = A + Bz + F 2 e iuJt - a2Z (15) 

Here the amplitude and decay lenght of temperature fluctuation in the sample 
depend in a nontrivial way on the diffusivity, heat conductivity of electron 
and phonon system and energy frequency of electron-phonon interaction. 
Equation (15) can be obtained from Eqs. (3) in the limit when P — > oo. In 
this situation the electron and phonon temperature are the same, the product 
P(T e — T p ) approaches to unknown finite value and the set of differential 
equations for P(T e — T p ) and T e = T p = T reduce to one, 

d 2 T ldT 

-d^ = am (16) 
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with 



where the solution is similar to Eq. (15) with the appropiate boundary 
condition at the surface of the semiconductor given by 

Q(z, t)\ z=0 = Q e (z, t) + Q p (z, t)\ z=0 = Q e + Q P + (AQ e + AQ P ) e lwt 

(18) 

At this point, the authors beleive that is important to remark the differ- 
ence between thermal waves and electromagnetic waves (for example). In the 
absence of rigorous theoretical guidance from first principles, several workers 
find it necessary to introduce arbitrary algebraic factor into their calculations 
in oder to get desirable fit to the data. As a consequence the controversial 
analogy between termal waves and electromagnetic waves arises. 

It is our contention that the main differences are the following: 

i) As pointed out by Mandelis 12 electromagnetic waves satisfy a hyperbolic 
differential equation while in the case for thermal waves satisfy a parabaloic 
differential equation. In other words heat conduction in solid is a diffusive 
process, quite different to electromagnetic wave propagation. 

ii) It is well known from standard electromagnetic books that the tangen- 
tial component of the electric field at the interface between two different 
media must to be continuous. In addition there are a reflected and trans- 
mited electromagnetic waves through the interface. However for the thermal 
waves there is not such similarity i.e., the electron or phonon temperature 
distributions in solid is not continuous in general at the interface of two layer 
system (see Ref. 10) only heat flux (energy conservation) is continuous at the 
boundary between two thermal different media. 

iii) It is also well known that a damped electromagnetic wave in a medium 
is of the form e *M- fe:E )- A;c where in general the decaying length staisfies the 
following inequality A" 1 > k^ 1 where k is the wave vector associated with the 
electric field, as mentioned before, information about electron and phonon 
thermal parameters can be obtained from photothermal experiments. In par- 
ticular when the electron-phonon interaction is strong such that k 2 e cu/a e , 
in this limit the decay length (3?<ti) _1 is smaller than the wave length (Sai) -1 
of the thermal wave associated with the electron system. This behavoir is 
impossible to get from the propagation of electromagnetic waves. This three 
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points show that the analogy between thermal and electromagnetic waves is 
not correct. 

4 Conclusions 

A theoretical analysis of thermal waves in thin films semiconductors has been 
studied. Using the appropriate boundary conditions, we obtain the electron 
and phonon temperature distribution taking into account the physical solu- 
tion for the static and dynamical parts of the heat-diffusion equation. For 
typical parameters of the heat conductivity el electrons and phonon in semi- 
conductors (k e ^> k p ), it is possible to obtain information about the phys- 
ical parameteres describing the diffusivity and electron-phonon interaction 
in photoacoustic experiments. It is shown, Eq. (12), that the photothermal 
signal is ultimately governed by the electron or phonon system depending on 
the relationship of the amplitude and the decay lenghts of the temperature 
fluctuation in each system and the sample size. We also showed that it is 
possible that decaying parameter of the thermal wave be smaller than the 
wavelength. However, it is well known that from Maxwell equations this 
situation is forbiden for electromagnetic waves. 

We have also derived exact solutions for electron and phonon temper- 
ature in semiconductors in the limit of weak and strong electron-phonon 
interactions. The above findings tell us that thermal waves in semiconduc- 
tors can propagate independent each other or propagate as a thermal wave 
in the sample with an effective diffusivity incluiding the contribution of the 
thermal parameters of the electron and phonon system. 
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